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THEORETICAL STUDY OF AIR FORCES ON AN OSCILLATING OR STEADY THIN WING 

IN A SUPERSONIC MAIN STREAM 

By I. E. Garrick and S. I. Rubinow 


SUMMARY 

A theoretical study, based on the linearized equations of 
motion for small disturbances, is made of the air forces on 
livings of general plan forms moving forward at a constant 
supersonic speed. The boundary problem is set up for both 
the harmonically oscillating and the steady conditions. Two 
types of boundary conditions are distinguished, which are 
designated “purely supersonic ” and “mixed supersonic.” The 
purely supersonic case involves independence of action of the 
upper and lower surfaces of the airfoil and the present analysis 
is mainly concerned with this case. A discussion is first 
given of the fundamental or elementary solution corresponding 
to a moving source. The solutions for the velocity potential 
are then synthesized by means of integration of the fundamental 
solution for the moving source. The method is illustrated by 
applications to a number of examples for both the steady and 
the oscillating cases and for various plan forms, including 
swept wings and rectangular and triangular plan forms. 
The special results of a number of authors are shown to be 
included in the analysis. 

INTRODUCTION 

This paper constitutes a theoretical study of the aero- 
dynamic forces on an oscillating or steady wing of finite 
span moving forward at a uniform supersonic speed. The 
treatment is based on the linearized theory obtained by 
considering only small disturbances in an ideal fluid. The 
wing is therefore considered to be a nearly flat thin surface 
at a small angle of attack and the flow is considered non- 
viseous and free of strong shocks. The theory in this case is 
equivalent to finding certain solutions of the wave equation in 
three dimensions with respect to a moving coordinate system. 

For the case of steady motion there exist a number of 
interesting solutions and methods. Among these may be 
mentioned the Yon Karman and Moore linearized treat- 
ment of slender bodies of revolution (reference 1), the Prandtl 
acceleration-potential method employed by Schliehting 
(references 2 and 3), the Busemann method of “linearized 
conical flows” (reference 4), studies of Jones, Puckett, 
Stewart, Brown, and Gurevich (references 5 to 9); and a 
method of Yon Karmen employing Fourier integral solutions 


of the two-dimensional wave equation and described by 
him as “acoustic oscillator method” (Wright Brothers 
Memorial lecture, Dec. 17, 1946). 

The corresponding unsteady or nonstationary problem for 
two-dimensional flow (infinite aspect ratio) may be con- 
sidered to be solved. In this connection there may be 
mentioned the work of Possio, Yon Borbely, Temple and 
Jahn, and the present authors (references 10 to 13). Of 
interest also are two wartime German papers by Schwarz 
and Honl (references 14 and 15). The corresponding steady 
plane case to which the nonstationary problem may be 
reduced is that treated by Ackeret. 

Results for the nonstationary or oscillating case are of 
great interest in the investigation of aircraft instability. 
The two-dimensional results have been applied to a study of 
flutter at supersonic speeds in references 12 and 13. Of 
more direct interest for this application are the three- 
dimensional results, especially for wings of swept plan form. 

The method used in the present study is to build up solu- 
tions of the equation satisfied by the velocity potential by 
superposition of the fundamental wave-potential solution 
for a spherical source. These solutions are also made to 
satisfy certain required boundary conditions on the airfoil 
surface. In the two-dimensional supersonic nonstationary 
case, which appears herein as a special limiting case, it can 
be proved that the procedure leads to a solution that is the 
unique solution of the given boundary problem. (For the 
problem of subsonic flow past a thin wing, reference may be 
made to the general treatment and method of Kiissner 
(reference 16) which also involves solutions of the wave 
equation.) 

Some qualitative features of the nature of the boundary 
problem may be mentioned here. Further remarks may be 
found in reference 17 and in Yon Karman J s Wright Brothers 
Memorial lecture. In the case of subsonic flow past an air- 
foil the whole field is influenced by the body. The concept 
of circulation has proved to be very useful and the Kutta 
condition has been used to specify the circulation by requir- 
ing smooth flow leaving the trailing edge. Thus, a deflecte_d 
aileron in subsonic flow influences the flow pattern over the 
whole wing even more importantly than over the aileron itself 
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In the case of supersonic flow the influence of the body is 
limited to only certain parts of the field of flow and generally 
the wake does not influence the upstream flow. The bound- 
ary problem for a three-dimensional surface moving at a 
supersonic speed can be classified into two' types referred to 
herein as “purely supersonic” and “mixed supersonic.” 
The definition of these terms is given in the analysis accord- 
ing to the parts of the field influenced by the airfoil, The 
purely supersonic case involving independence of action of 
the top and bottom surfaces and no reflecting surfaces in the 
field. Thus, in the purely supersonic case, a deflection of 
the aileron would produce only a local effect near the aileron; 
in the mixed supersonic case, it may have a decided influence 
on the part of the wing adjacent to the aileron or on other 
parts of the wing. For a given wing both types of problems 
may be involved. 

The treatment used for the purely supersonic cases, involv- 
ing source and sink distributions to account for the action 
of the body, is believed to be exact within the framework of 
the linearized theory. The upper and lower surfaces of the 
airfoil are regarded as acting independently, each surface 
being “unaware” of the presence of the other. The treat- 
ment is thus analogous to that of sound in a moving medium 
generated by the motion of pistons imbedded in an infinite 
plane. This flow picture is obviously incomplete in the 
mixed case and more complicated distributions (doublets) 
are also required. For some purposes, however, the simpler 
treatment may still be used in conjunction with appropriate 
correction factors. Also, for steady flow past a symmetrical 
airfoil at zero lift, the simpler treatment can be employed 
for study of the wave drag. 

The object of the present paper is to develop the expres- 
sion for the velocity potential in the purely supersonic case, 
based on the elementary solution for the sound source 
moving uniformly at a supersonic speed, and to indicate 
its application by a number of special examples. 

SYMBOLS 

<j> disturbance-velocity potential 

x', y', z’ rectangular coordinates for fixed system 

x, y, z rectangular coordinates attached to source mov- 
ing in negative ^-direction; also represents 
field point being influenced 
rectangular coordinates used to represent space 
coordinates in fixed system 

£, 17 , i rectangular coordinates used to represent space 

location at source distribution A ( |, 17, f) 
t, T, t' time 

v velocity of main stream 

c velocity of sound 

M Mach number (v/c) 

r distance defined by equation (8) 

n, r 2 time function defined in equation (7a) 

g function defining airfoil surface ( y—g(x , z, t)) 


foj fi, 1 * 2 , fi limits defined in equation (10) 

9 variable used instead of f defined by relation 

preceding equation (15a) 
p pressure 

p 0 reference pressure 

p density 

a angle of attack 

a time derivative of a 

co angular frequency 

w(x,z,t) vertical velocity factored in equation (14) as a 
space function W(x,z) and time function w(i) 
A angle of sweep 

h vertical displacement 

li time derivative of h 

ANALYSIS 

WAVE EQUATION AND SOURCE SOLUTIONS 

In the linearized theory based on small disturbances the 
equation satisfied by the velocity potential for the propa- 
gation of sound waves of small amplitude is the wave equation 

1 bV b 2 <£ , 5 2 (£ , 

c 2 dt ' 2 bx'^by'^bz ' 2 ^ 

The fluid medium is considered at rest at infinity. 

In the treatment of linear partial differential equations the 
so-called elementary or fundamental solution is of great 
importance since general solutions can be built up by dis- 
tributions of elementary solutions. From a physical point of 
view the elementary solution may correspond to a source. A 
discussion of the nature of element ary solutions for hyperbolic 
differential equations of a general type has been given by 
Hadamard (reference 18), who makes the cardinal statement 
that “every result of the theory can bo and has to be deduced 
from the consideration of the elementary solution only.” 

A fundamental solution of equation (1) from which general 
solutions may be formed is that of a source of sound fixed 
in the medium 

*0 (2) 

where 

In equation (2) the fixed source is located at the point 
(!', n', f0, the strength of the source is A(£', y', I') fit'), 
and the minus sign indicates that the spherical waves are 
diverging from the center of the disturbance. 

Another closely related solution of equation (1) is that of 
a fixed point source for which the spherical waves arc con- 
verging onto the source 

(3) 

The wave potential in equation (2) is often designated 
“retarded” and that in equation (3), “advanced.” 
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It is intended to consider thin lifting surfaces of small 
curvature which are moving forward at a constant supersonic 
velocity v and which may be performing small oscillations 
normal to the direction of v. The direction of v will be that 
of the negative r-axis and the surface will be replaced by a 
distribution of moving sources in the rs-plane (fig. 1). 



Figure I.— The coordinate system used. Thin lifting surface in xz-plane moving at a 
constant supersonic velocity v in the negative ^-direction. 


Consider a source moving in the negative r-direetion with 
uniform velocity v and a rectangular coordinate system 
attached to the moving source. If the new coordinates are 
designated by x, y, z, t, where x—x'-ri't', y=y', z—2', t=t', 
the equation satisfied by the potential is 


or 


I /i> 1 , AY ■ 

dr/® dr 2 ‘ by^bz 2 

J_ ^ 4. A d 2 d , /V_ \ d 2 d b"-<p b 2 <b 
c 2 bt 2 ^ c l bxbt 1 \c 2 / dr 2 by 2 bz 2 


( 4 ) 


This equation is satisfied by the potential of sources of sound 
in motion through the medium with uniform velocity v in 
the negative r-direction. It is also the equation satisfied by 
the disturbance velocity potential for a fixed body creating 
small perturbations from an oncoming main stream of 
velocity v in the ar-direction. A brief derivation from_ hydro- 
dynamical principles is given in appendix A. 

It is known from the classical study of the wave equation 
(reference 16 ) and can he verified by direct substitution that 
a solution of equation (1) is transformed to a solution of 
equation ( 4 ) by means of the following substitutions, corre- 
sponding to a combination of the Lorentz transformation 
and a Galilean transformation: 


VI' -M* 

y'=y 

z'=z 




Cyl—M 2 I 


( 5 ) 


where M, the Mach number of the main flow, is v/c. 

For the purpose of studying the supersonic case (M)> 1), 
it is more convenient to employ modifications of transforma- 


tion ( 5 ) obtained by multiplying the right-hand side by tbe 
constant 1 fij 1 — M 2 , or 


x = 


x 

1— M 2 
V 

TOOT 2 


V 

t' = t-L 


/l-M 2 


xM 

r c(l-M 2 ) 


( 5 a) 


The particular solution of equation ( 4 ) that corresponds 
to a moving source will be seen in the following discussion 
to be analogous to a solution of equation (1) given by the 
sum of potentials in equations ( 2 ) and ( 3 ), namely, to 


<Po- 




( 6 ) 


The desired solution of equation ( 4 ) corresponding to equa- 
tion (6) is obtained with the aid of the substitutions ( 5 a) as 


, _AV ( M r\ , / M x-Z 

^~r\jV cAP-1 c) +J V c M 2 - 1 


( 7 ) 


where 


1 

~AP— 1 


■J{x-W-{M 2 ~l)[{y-nY-+{z-m ( 8 ) 


(The term -yjl—M 2 in equations ( 5 a) causes no difficulty 
since only the squares of the space coordinates are needed.) 
This solution for 4 >o may be expressed in the form 


<Po— y L/(* — v 2 ) -h A (f — ri)l ( 7 a) 

where 

M x—£ , r 

T2 c M 2 -l + c 

_ M x-l r 
11 c M 2 — 1 c 


and where r is defined as in equation (8). The constant 
A(£,i?,r) could of course have been included in the functional 
symbol / but has been separated for convenience. It may 
be considered to represent the space variation of the source 
strength as distinguished from the time variation of strength. 
For a moving source of constant strength the time function 
may be considered equal to unity and the potential expressed 
as (reference 2) : 



It will be recognized that the solution, equation ( 7 a), is 
valid in a conical region, the so-called “Mach cone,” opening 
aft of the moving source. Outside of this conical region, 
defined by the equation r=0, the flow is undisturbed. 
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The result expressed by equation (7) may be considered 
physically from two points of view. In one, as considered 
by Prandtl (reference 2), a source of variable strength 
moving along a certain path is replaced by a continuous 
succession of fixed source pulses distributed along its path 
acting consecutively one after the other. Each pulse, con- 
sidered fixed in an absolute coordinate system, emits a 
spherical wave traveling at sound speed and the coordinates 
of the center of the spherical surface are £+»£, »?, f. The 
radius , vector R of a point (x, y, z ) with respect to this 
center is 

R= Vfr- (*+rf)] 2 + (y-n) 2 + O-f) 2 ' 

The time at which the spherical wave passes the point (x, y, z) 
is 



Eliminating R between the preceding two relations results in 

c 2 t 2 - (x-l—vt) 2 - ( y—v ) 2 — (z— 1) 2 =0 

The roots of this quadratic equation in t are precisely the 
quantities n and t 2 defined in equation (7a); that is, the 
field point (x, y, z) is influenced at time t by two waves which 
originated at times r 2 and it earlier. It is of interest to 
observe that, in the supersonic flow, both roots are real and 
positive and have physical significance; whereas, in the 
subsonic flow, only one root is positive and of physical sig- 
nificance, In the supersonic case the field of influence of a 
source is the particular Mach cone with vertex at the source, 
■ and through each point in this region at instant t, there pass 
two spherical surfaces representing the waves originating at 
times n and r 2 earlier (fig. 2). 



Figure 2.— Field, of influence of spherical source moving at a constant supersonic velocity. 
Point (a:, jf, z) is affected at time t by two pulses originating at times n and ti earlier. 


From the other point of view of the result (equation (7)), 
a single diverging spherical wave-pulse is considered. Let 
this wave originate at the point (f, q, f) at a time T (fig.” 3) 
and consider its effect at a point (x, y, z ) (within the Mach 
cone whose vertex is at point (£, ij, £)) moving with a velocity 
greater than that of sound. Clearly at a later time T+n 
the moving point penetrates the wave front and at a still 
later time T-tt 2 it emerges from the wave front. The 
potential at point (x, y, z) changes only on entering and on 
leaving the wave front and the two terms in equation (7) 



Figure 3.— Influence of single spherical wave pulse. Spherical wave originating at points 
(£, v, f) at time T influences point ( x , v, 2) fixed relative to ($, ij, f) and moving at a con- 
stant supersonic speed, at times n and n later. 

correspond to these two effects. The factor 2 appearing in 
the potential for a constant source moving at a supersonic 
speed also has its origin in this physical fact, in contrast to 
that for a source moving at a subsonic speed, where the field 
point penetrates the wave front hut never emerges and where 
the corresponding factor is unity. The two-dimensional 
supersonic case involves cylindrical waves and the potential 
of the . point (x, y) is continuously changing from the 
time the point enters to the time it emerges from tiro 
wave (reference 13). Observe the interesting geometrical 
property of r (equation (8)); namely 2 r is the difference of 
the radii of the spherical wave at time r 2 and at time n, that 

is, (t 2 — it). (Observe also that the potential which for- 

mally appears in equation (6) as the sum of potentials, half- 
advanced and half-retarded, transforms in the moving 
coordinates to a sum of retarded potentials in which the 
original retarded part is associated with the diverging 
spherical concave wave from which the point is emerging 
and the original advanced part is associated with a diverging 
convex wave into which the point is penetrating.) Recent 
papers of interest in connection with moving acoustical 
sources are references 19 and 20. 

SURFACE DISTRIBUTION OF SOURCES 

Sources and sinks of the type 4> 0 will now be distributed to 
represenfijhe upper and lower surfaces of a thin airfoil. 
The procedure to be followed is that used in the two- 
dimensional case (references 10 to 13) where the upper and 
lower surfaces are considered separately. Also the total effect 
may be separated into an effect of the mean-camber surface 
and an additive effect due to thickness alone. In most of 
the applications, unless stated to the contrary, the mean- 
camber surface is considered. 

Let a. continuous distribution of sources be given over the 
mean-camber surface. The airfoil is considered so thin and 
flat that the source distribution may be treated in the 
arz-plane (fig. 1). The airfoil surface may be considered 
moving at a constant speed v in the negative z-direction (or 
fixed in a stream moving in the ^-direction). The effect at 
a point (a;, y, z) at time t of a distribution of sources of posi- 
tion magnitude A((-, 0, f) is given by an appropriate 
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integration, over a region of the -plane of the form 


4 ,{x,y,z ! t)=ff 4 > 6 dtd£ (9) 

where <po represents the function given in equation (7) with 
17=0. 

The total effect at the point (2, y, z) is the sum of the 
effects of all disturbances having their origin within the 
Mach cone with vertex at point (2, y, z) and opening in the 
upstream direction. This conical region need not extend 
into the undisturbed part of the flow; that is, it need not 
extend beyond the most forward surface envelope of the 
Mach cones of influence of the body. There are essentially 
two types of boundary conditions that need to be dis- 
tinguished, designated by the terms “purely supersonic’ 1 
and “mixed supersonic.” A point of the boundary belongs 
to a purely supersonic case if the upstream facing Mach cone 
contains, in the part of the 22-plane not considered occupied 
by the body, no disturbed fluid having a component normal 
to this surface. Otherwise, the point belongs to the mixed 
supersonic case. A sufficient (but not necessary) criterion 
for the purely supersonic case is that the component of the 
main stream normal to any edge or contour of the plan form 
in the 22-plane (contained within the upstream facing Mach 
cone of the given point) shall be supersonic. There is no 
downwash ahead of the body, no holes are in the body, no 
spilling of fluid occurs around edges, and no reflecting surfaces 
are in the flow field. In this case the upper and lower sur- 
faces of the airfoil are considered to act independently of 
each other; a disturbance created on one side does not affect 
the opposite side. The flow can be considered to arise from 
the appropriate movement of small pistons acting at the 
regulating or generating surface. This condition is in con- 
trast to that of the mixed supersonic case, for which the 
effect of the disturbance spills over the edges or sides, and 
a disturbed fluid region (downwash) may exist ahead of the 
body. Thus, points of a triangular surface, moving vertex 
foremost and completely outside of the Mach cone associated 
with the vertex, belong to the purely supersonic case. If 
the triangular surface is inside the Mach cone associated 
with the vertex, the points belong to the mixed supersonic 
case. Of course, for a given surface, both cases may be 
involved. A few examples are shown in figure 4. 

In the purely supersonic case the circulation concept plays 
no particular role and the drag associated with lift or thick- 
ness may properly be denoted as wave drag. In the mixed 
case the flow retains subsonic features and the drag associated 
with the lift is sometimes denoted as induced drag. 

Although the treatment given for the purely supersonic 
case is believed exact within the limitations of the linearized 
theory, an exact treatment of the mixed supersonic case is 
not available. These problems involve greater difficulties 
in the boundary conditions, for the flow to a certain extent 
acquires features of a subsonic flow in that the fluid field 
“senses” the approach of the body. Thus, in certain cases, 
conditions at the leading edge, at the trailing edge, and in the 
wake must be specially taken into account. For some pur- 
poses and in certain problems, however, it may be useful to 
treat the mixed supersonic case in the same manner as the 





i 



Figure 4.— Illustrations of plan forms with “miied supersonic” regions (unshaded) and 
] ; “purely supersonic” regions (shaded). Mach lines are shown dashed. 


purely supersonic case and to introduce appropriate correc- 
tion factors. 

The region of integration in equation (9) is the part of the 
body (in the £f-plane) cut out by the upstream opening 
Mach cone with vertex at point (2, y, z). This region in 
general depends on the plan form of the body as well as on 
the point (2, y, z). With the understanding that the leading 
point of the body is at £=0, the integration may be written 


where 


<p(x,y, s, t) = 4 >o dg d£ 


ri=z-r<> 


r ° yM 2 -l y 
^=x—y^M 2 —l 


( 10 ) 


The limits of integration and f 2 in equation (10) may be 
recognized as the distances from the £-axis to the near and 
far sides, respectively, of the hyperbola defined by the inter- 
section of the cone r=0 and the plane 17 = 0. Thus, from 
equation (8), with ij=0, fi and are recognized as the roots 
of the equation 

M,M)=^=3y V(f-iTi)(f*-r) = 0 (11) 


The limit in equation (10) represents the ^-coordinate of 
the vertex of the hyperbola and is defined by the condition 
f l= f 2) that is, by l"o=0. The point (fi, z) is the farthermost 
downstream point which can affect the point (2, y, z) 
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BOUNDARY CONDITION 


The strength of the distribution of singularities in equa- 
tion (10) will now be determined by the boundary condition 
of tangential flow along the airfoil surface. The boundary 
condition may be expressed as 


&4> 

by 


y =0 


=w(x, z, t ) 


-_ 0 

0 bx^bt 


( 12 ) 


where the airfoil shape is defined by y=g{x, z, t) and where 
the two terms represent the normal velocity induced by the 
airfoil shape and by its own proper motion. It is shown in 
appendix B and can also be made clear by physical reasoning 
that as y approaches zero from the positive side (y— »+0) 


bcj> 

by 


= -2x(M 2 -l)A(x,0 ) 2)/(f) 


or, briefly, 


A(x,z)j{t) = - 


' 2ir(AP— 1) 


w(x, z, t) 


(13) 


As y approaches zero in the negative half plane, an equal 
and opposite result is obtained. Equal source distributions 
on the upper and lower surfaces therefore result in a discon- 
tinuous vertical-velocity distribution . near the., plane y = 0 
and may be used to represent symmetrical thickness distri- 
butions. The source distribution representing a thin body 
with arbitrary thrclaress"disti'ibution is in general unequal 
on the two surfaces. The effect of thickness es discussed in. 
a separate section. A representation of the mean-camber 
surface alone may be obtained by placing equal and opposite 
sources on the under surface in. proximity to the sources on 
the upper surface. The potential <j> is to be understood in 
the subsequent analysis to be prefixed by a ± sign, plus for 
the upper surface and minus for the lower surface. The 
vertical velocity will in general be measured positive upward. 

It is convenient to express the vertical velocity in equa- 
tion (13) in separated form 


where 


w(x,z,t) = W(x,z)w(t) (14) 

W(x,z) = — 2tt(M 2 - \)A(x,z) 

SURFACE POTENTIAL 

The total potential for y = 0 may now be expressed by 
means of equations (10). and (14) as 


cx rti 

4(x,z,t)= \ bo d’y dk 
Jo JSi 


. w(t— rj+wjt— r 2 ) 


■5&C r w^m=t ds d( <i5) 


where, for y= 0 (see equations (7a), (10), and (11)), 


V(t-r,)(r 2 -r) 

c/3 2 cf 3 


Ur 


M(x-Q V(f-ft)QW) 
cf. 3* 0 


to 

(* 2 =2+fo 
.2— I 


To— - 


J3 


j9=VH2*-l 

and where it is understood that U r (£, f) = 0 at any point off 
the body or where the integrand is not real. 

Equation (15) may be put into a simpler form by substi- 
tution of a new variable 9 instead of f, which is obtained 
from the relation (see appendix B) 


or 


2 f=(f 2 — fi) cos 0 +f 2 +Ti 
f=fo cos 9+z 

The surface potential (equation (15)) may thenbe written as 

<f>{x, z, t) — ~ 2 J o IE[h f(5)][w(i— Ti)+w(t— r a )] d.8dt 

(15a) 

where 

Ti= X -0 (M- sin 5) 
r 2 =^ (M+sin 6) 


y- 

=cos _1 — 


x—i 


p 


Equation (15) represents the central result of .the analysis 
and within the limitations already discussed may bo applied 
to wings of any plan form in steady motion or performing 
small oscillations. In the stationary or steady case, <t> docs 
not depend on time and the function w(t) is to be replaced 
by unity. Then, in equation (15), w(t— n) r 2 ) is to 

be replaced by 2. 


PRESSURE RELATIONS 


For the sake of reference, relations for the pressure and 
the lift and drag forces are given here. The disturbance 
pressure (local static pressure minus the pressure in the 
undisturbed stream) may be written as 


v=- 


d<j> 
p dt 




(10) 
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The pressure difference (positive if acting downward) at any 
point (x, z) may be expressed as 


The integration with respect to 8 may be readily performed 
with the aid of the relation 


AP=Pu—Pl 

where the subscripts U and L refer to the upper and lower 
surfaces. For the mean-camber surface Pl=—Pu and 

AP=-2p(^+»§i) (17) 


The total forces on the airfoil in the y-direction and 
^-direction are given by 

F=Lift= S Sv dx dz 

A=Drag = —ffp dy dz 


where the integration is to be taken over the complete airfoil 
surface. Expressed as integrations over the plan form 


7= 

X= 


I'Jc p L ~Pu)dxdz 

/Mix-®: 


(18) 


\dxdz 


Finally 


where 


— f cos (\ sin 9)d6= J 0 (A) 
* Jo 


6(x,t) = — 




Ut,x)=e 


x—j M 

c * Jo 



( 20 ) 

( 21 ) 


This result for the velocity potential is identical with equa- 
tion (11) of reference 13 and is used therein as .a basis for 
calculation of the nonstationary two-dimensional case. 

In the steady case, o=0 and /(?, x)=l. The expression 
for the velocity potential is 



( 22 ) 


where 


ir(0=® 


dy. 

dx 


This formula or the pressure relation 




i? dy 
bx~ p (S dx 


It is often convenient to separate the slope terms as follows: 






= cc+<r L 


where a is the conventional direction of the main stream with 
respect to a reference chord, and tr u and <r L are the local slopes 
of the airfoil surfaces measured with respect to the reference 
chord and positive in the same sense as a. 


APPLICATIONS 

WING OF INFINITE SPAN AND ZERO SWEEP 

For the first application of equation (15) the results for 
both the oscillating and steady two-dimensional case will be 
derived. For the harmonically oscillating wing having 
identical motion in every chordwise section, the vertical 
velocity can be written in the complex form 


Then 


w(x, t) = W(x)e iat 
w(t— r t ) +w(t— r 2 ) =e iu ‘ 


=e iut ( 2e 


Equation (15a) becomes 


2e 2 cos 6>- r2 ~ ri ^ 


*(*, t) = ~~ J q J o cos j 2 sin 8) ddd^ 

^ ' (19) 

where (3= 1. 


applied to both the upper and lower surfaces of the airfoil 
leads to all the results of the Ackeret theory. 

WING OF INFINITE SPAN WITH ANGLE OF SWEEP 

Consider an infinite wing with angle of sweep A (fig. 5), and 
assume that all sections in the flight direction are identical 
in shape and that the wing is undergoing harmonic motion. 
In general, the vertical velocity is can be written in the com- 
plex form 

w(x, z, t) — W ( x , z)e iul 



Figure 5. — Sketch, for ’sing of infinite span with angle of sweep showing region of 
Integration (shaded). 
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If each section normal to the leading edge is performing the 
same motion, the form of W(x, z) is W(x—z tan A). If the 
wing is assumed to perform pure vertical motion alone, then 
W{ x, z) is a constant. If the wing is assumed to rotate about 
an axis x= Constant, then W(x, z) is of the form W{x). 


The potential is of the form (fig. 5) 

4>(%, M) = Jir Fde d£— £/; Fde d£ 

where 

F (£, e, t ) = W& f (fl)) [e^'-rj+e*^] 

and where 


, a:— g ff 

1-/3 cot A 


(23) 


_ X+Zf 3 

^ 1 + /3 cot A 

/£ cot A — Z -\ 

9, = C0S ( x _ ( ft) 

The values of the limits £ 3 and £ 4 are found by solving for £ 
in the relations fi=£ cot A and f 2 =£ cot A which represent 
the intersections of the Mach lines through x with the lead- 
ing edge. The limit 6=d } corresponds to £=£ cot A, the 
leading-edge line. 

When IF(£, z) is a constant or a function of £ only, the 
velocity potential can be expressed as 


4>(X, 2, = [Jj* W r (£)J(£, X, ,)*] 

(24) 

where /(£, sc) is as defined previously and 

/i(£, z)=\ p cos(^ ^ sin 6^ dd (25) 


Observe that the integral involved in equation (25) for i) 
reduces to the Bessel function of zero order when 8i = tt as 
in equation (19). This interesting integral may therefore 
be called an “incomplete” Bessel function of zero order. 
Systematic investigation of its properties would appear to 
be desirable. 

For the infinite swept wing in the steady case the frequency 
c<j may be made equal to zero in equation (23). Consider 
as a simple example the case of a thin wing at a small constant 

angle of attack a, that is, ^=— a. Let the angle of sweep 


be less than the complement of the Mach angle, that is, 
/3 cot A>1. (Otherwise the case involves the mixed- 
supersonic-flow conditions.) From equation (24) with co=0, 

/(£, x)=l, and / t (£, x, z)=% 

' 7T 


^ =7 [/I «-aV C c ° t 4 s ) * 


=Va 


x cot A—g 
V/3 2 cot 2 A— 1 


(26) 


The local pressure difference is given by 
^ _ 2 pv i a (3 cot A 

p ~ p TFc otHT-i 


( 27 ) 


This equation reduces, for A = 0, to the Aekeret result 

2 pv t a 


Let the index n refer to quantities measured normal to the 
leading edge. Then 


and 


a n =a sec A 
v n =v cos A 

M n =—=M cos A 

_ _ c 

2 pVnOCn 


a result similar in form to the expression for p 0 and already 
stated by Busemann (reference 17) in 1935. (See also 
reference 6.) 

The harmonically oscillating case with W(x, z) assumed 
to be of the form W(x—z tan A) leads in a similar manner to 
a result analogous to equation (20). 


RECTANGULAR WING OF FINITE SPAN (ZERO SWEEP) 


Consider a harmonically oscillating rectangular wing of 
finite span as in figure 6. Region I is described as purely 
supersonic and region II as mixed supersonic. The higher 
the aspect ratio and the stream Mach number, the relatively 
smaller the region II becomes. 

The potential for region I for identical motion of each 
chordwise section is exactly that given for the infinite wing 
in equation (20); however, more general types of motion 
involving spanwisc variation may also be treated. For 
example, let the wing perform harmonic oscillations in 
vertical bending and in torsion about a spanwise axis' x=x Q 
in certain prescribed spanwise modes. Then, with a and 
h used to describe angle of attack and vertical position 
(%• 7), 


a=a 1 (z)a 2 (t) 1 
h=k l (z)h 2 (t) ) 


(28) 


where a, ( 2 ) and hi(z) represent spanwise modes and 

a 2 (t)—a a e iut 

h 2 (t)=h 0 e i “ t 


and a 0 and h Q are constant complex amplitudes. The vertical 
velocity (w measured positive upward, k positive down- 
ward) may be expressed as 


w(x, — (x—x 0 ) a] 


(29) 
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Figuee 7.— Sketch for oscillating rectangular wing (ft is positive downward, 
a is positive clockwise). 


With the use of equations (15a), these potentials may be 
expressed as 

it j: -j; c ° s (h sin e ) d& ^ j 

fo = Zpj 0 Ai(r)cos^si ndjdddz j (31) 

<Pi=|| r e _i3 £ (f— z 0 )«i(r) cos sin 8^ ddd£] 
where 

wilf(ar-€) 

s_ c(M 2 -l) 

and, expressed as a function of 6, 

ai(f)=ai(z+fo cos 8) 

£i(f)=^i(z+fo cos e) 

If the modal functions in equations (31) are a E =A, = 1, the 
potential corresponds to that giveD by equation (20) for the 
two-dimensional case. (See also equation (14) of reference 13.) 
It is of interest to consider modal functions for a £ and 
A l of the type (f/«)“ where s is the semispan. For modal 
functions of this form the typical integral involved in 
equations (31) may be expressed as 

F n = ( 2 + To cos 0) n cos (jj siEL 0 ) dd ( 32 ) 
With the substitution of 8 for 8, F„ may be written as 


F = 


J 2 (z+ fo sin 8) n cos cos ® 
J ^( z ~ fo si 11 ®)“ cos ^ijj-cos 


de+ 

de 


The further reduction of F is made with the aid of the fol- 
lowing relation (reference 21) : 



cos 8) sin 2i 8 dd 


For example, the case n=Q corresponds to constant modes 
and yields for the potential the result already given by 
equation (20). The case n=l corresponds to linear modes 
and the function Fi becomes 


Let the potential (equation (15)) be separated into the 
form 

4> = 4>a J r4>s J r < i>* ( 30 ) 

where the various <j>’s are associated with the corresponding 
variables in equation (29). 


This relation utilized in the equation for the potential yields 
a result that is the two-dimensional-case result multiplied 
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by the factor 2 /s. The case n=2 corresponds to parabolic 
modes and the function F 2 becomes 

F ! =.«(j0-W>-8.7,(£) 

When F 2 is used in equations (31), the J 0 term yields an 
integral of the tj r pe given by equation (20). With the use 
of the relation JJ\) = — J/(X), the Ji term also yields an 
integral of the same form. This type of reduction to the 
form of equation (20) maybe made in general for any integral 
index n by means of the recurrence formulas for Besselfunc- 
tions, and thus use may be made of the numerical procedures 
used for equation (30). (See reference 13.) 

It may be of interest to treat the potential for the mixed 
supersonic region II (fig. 6) as though.it were part of a purely 
supersonic region. The equations corresponding to equa- 
tions (31) are 

0" e ~io a ^ C ° S G? Sin 6 ) dd d *~ 

C0S (If sin 0 dd ^ <- 33 ) 

and similar equations for and 4>„. The limit | 2 is found 
as the value of | for which f 2 =s or 

& =x—(s—z)P 

The limit 0=0 corresponds to and the limit d—0 2 

corresponds to £= s or, from equation (15a), 

»,=cos-' (0^) 

The last term in equation (33) leads to integrals of the 
“incomplete” Bessel function type as mentioned for the case 
of the infinite wing with angle of sweep. 

The foregoing results for the oscillating rectangular wing 
will now be specialized to the steady case (u=0, g—0, 
ai (f)=l, a 2 (t)=a, the constant angle of attack). Then, 
from equations (31), the velocity potential for region I, is 

<t>a=JX (34) 

For region II, from equation (33), 


It may be observed that p is constant along rays from the 
£ 2 

tip =Constant. Along the ray corresponding to the 

g 2# 

Mach line from the tip, 1 and p takes on the constant 

value p x . Along the ray corresponding to the tip z=s, half 
of this value is obtained. This edge condition is physically 
incorrect since the assumption of the independence of the 
two surfaces of the airfoil is not correct near the tip. 

This particular problem has been treated by Busomann 
(reference 4) by his method of conical or perspective sym- 
metry. The condition along the ray corresponding to the 
tip is p = 0 and Busemann’s result for region II is 

2 pv-a 1 /, 2(s — z) „Y 

l 1 — —») 

The total lift over region II is one-half of that of an equal 
area of region I. A comparison of this result and equa- 
tions (36) is shown in figure 8. This comparison gives an 
indication of the errors involved in the assumption of inde- 
pendence of the two surfaces near the rectangular tip and, 
conversely, it gives an indication of the appropriate correc- 
tion factors required to allow for the tip effect. It appears 
that equations (36) overestimate the lift over all of region II 
2 

by a factor I—- or by approximately 36 percent. Actually, 

equations (36) apply to the edge of a rectangular wing 
adjacent to a straight surface barrier at zero angle of attack. 

THICKNESS DISTRIBUTION 

It has already been remarked that the treatment em- 
ployed in the analysis mainly for the mean-camber surface 
can also be applied to obtain the effect of thickness. In 
equation (15) the vertical velocity IF(f, f) may be specified 
for both the upper and the knver surface. 

As an example, consider a plan form such as that shown in 
figure 9 in steady supersonic flow. Let the airfoil section 
shape, for convenience chosen symmetrical and independent' 
of span, be defined in the center section by y=g(x ) (and in 
any other section by y=g{x — 2 tan A)). Then, for the upper 
surface, 

W(x, z)=vct-\-vg' 
and, for the lower surface, 




va 

T 



The actual integration in equation (35) may be easily per- 
formed but is not required for the purpose of obtaining the 
local pressure. 

The local-pressure difference is directly obtained for 
regions I and II from equations (34.) and (35) as 


Pi 


2pv 2 a 

& 


Pn~ 


2 pv 2 a 
/3 



(36) 


W{x, z)=va—vg ' 

where g' is the derivative of g with respect to its argument. 

The velocity potentials in the various regions in figure 9 
are of the form 

n * rzi re 1 v 

— ax- d * + ££ Fdd rf£ ) * 

*m=*i-(J 6 ° J ft Fde Fd8 


( 37 ) 
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Figure 8.— Pressure for rectangular edge as given by Busemann and by the approximation 
given in equations (36). 



Figl'ee 9. — Sketch for swept wing for the case ,4>i showing regions I, II, III. 


f 2 =| cot A; and | 5 is the value of | for the leading edge of 
the tip. Then 

, i=c03 -.(1MA=£,) 

d 2 = cos -1 - — f j3 
x — | 

. x—zfi 


& 


2+2$ 
1 + A 


$6 


X—Z& 

1 + A 


la =x—(s—z)p 
| S =s tan A 
A=p cot A 

If, for example, the distribution function F is a constant K, 


<i>i=K 


X cot A — 2 
A 2 -l 


k r *« k r«‘ 

*u=fc-- Ml+- 0& 

n 5 n Jts 

k r° r** 

<t>ni=<Pi--\ (t- 0,)<?$ + (x-0 4 )<?| 

7 J& Jo 


The corresponding local pressures are 
P vK A 


Pi=- 


P V A 2 — 1 


pi’-ST ki / 1 , l-f-J-BX 

3, “~- r3JCit 1 W cos, XFB-; 


Pm — " 


J 1/ ^1+AC 

P yj 2 - iA 


cos ~‘XfC +cos A 


1 -AG\ 
A-C ) 


ptK A 2 , 


i /4!z 

Vi- 


— 1 

<? 2 


where 


A= p cot A 
(s—z)p 


B 


x—s tan A 


c=-/s 

x 


The limits in the foregoing integrals are as follows: the 
limits 8 = 8 l and 5=5 4 correspond, respectively, to the leading- 
edge lines f=| cot A and f= — | cot A; 5=0 and 5 =tt 
correspond, respectively, to the Mach lines f=f 2 and 
i = Sii la and fs are obtained, respectively, from the relations 
fi = f cot A and f t = — | cot A; 5 = 5 2 corresponds to f=s; | 2 is 
obtained from the relation f 2 =s; li is obtained from 


The constant K may be interpreted as ca associated with 
constant angle of attack. In this case, region II is to be 
regarded as a mixed supersonic region and the result given 
is not the appropriate solution for this region. If the con- 
stant K is interpreted as g'= Constant, the results are 
applicable to a thin symmetrical wedge of half vertex angle 
K and may be employed to yield the wave drag according to 
the linearized treatment. 
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Jones (reference 5) treats symmetrical airfoils of various 
plan for ms at zero lift use of pressure potential. The 
use of velocity potential leads to the same results as given 
in reference 5. Thus, equations (13) and (14) of reference 5 
for a wedge correspond to the preceding results. The 
velocity potential in general is more useful to treat pressure 
distributions for a given body; whereas, the pressure poten- 
tial may be more readily adapted to treat airfoil shapes and 
plan forms associated with desired types of distributions 
of pressure. 

TRIANGULAR PLAN FORM 

The triangular wing (fig. 10) extending across the Mach 
lines from the vertex may serve as a final example. For the 




Figure 10,— 1 Triangular wing in a supersonic stream and pressure distribution. Case 
sketched corresponds to A=3G°, 


steady case of a vanishingly thin surface at angle of attack a, 
the velocity potentials and pressure relations for regions I 
and III are equivalent to those just discussed in the preceding 
section. The lift AL on a strip A* of the triangle located at 
abscissa x from the vertex is given by 

Cx cot A 

AL= Ax A p dz 

J —x cot A 


4pv-a. 

IA-1 

4pV 1 a. / 


L (3 2 A- \ 

1 A+l 

h /3 s A V 

1 Ya+i/J 


where The two terms correspond to the integrations over 
regions I and III, respectively, and where A— l 3 cot A. Then, 

. T 4pv 2 a , . 

AL= M 2 - 1 A x Ax 


The area of the strip is 2x Ax cot A, and hence the lift coeffi- 
cient is independent of x and equal to 

ri AL . 

L ~fl M2Ax Ax\ 

4a 

~T 

Gurevich (reference 9) treats this case, and his relations can 
be shown to be equivalent to the foregoing ones. The pres- 
sure distribution is illustrated in figure 10, where po, the 
reference pressure, is 2py 2 a/j3. Observe that the pressure 
area above the unit ordinate cancels the area of pressure 
deficiency below the unit ordinate. Also shown in figure 10 
is the distribution of pressure as the half vertex angle of the 
triangle approaches the Mach angle. 

The triangular wing inside the Mach cone from the vertex 
requires _a more elaborate treatment (references 7 to 9). 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., June 4, 1947. 



APPENDIX A 


DIFFERENTIAL EQUATION FOR THE VELOCITY POTENTIAL 

A derivation of equation (4) is given briefly here. The 
condition for irrotational flow is 


curl v =0 


(Al) 


and this relation implies that a scalar velocity potential <t> 
exists, such that 

v=grad 4> (A 2 ) 


The general equation of continuity 
Sp_ 

may be written as 


bt 


+div pv =0 




(A3) 


where differentiation following the particle is denoted by 
g=|+(v-grad) 

and V 2 =div grad is the Laplacian operator. 

From Euler ^"equations, or from the general Bernoulli 
relation, 

^ r ..2 r* 

(A4) 


^ + L 2 + f 

& 1 2 \) p U 


where a spaee~constant function of time has been included 
in <f>, and wheredt has been assumed that p is a function of 
p“onlv. With the use of equation (A4) and the acoustic 
relation, 

c ~~dp 

where c is the local variable speed of sound, it follows that 

grad (^f + 0 = ~J gl ' ad ^ 


= — - grad p 


and 


/b<f> _ 1 _ dp 

bt\bt + 2 )~ p bt 


£2 dp 
p bt 


With the aid of these two relations the first term in equa- 
tion (A3) becomes 

., 2 \ 

i / U W . UV . T I 

. i * . I I _ t r W 

p Dt 


l Dp 1 /b 2 d> Sr 2 , v 2 \ 

-tf— F+ 5T+ 7 -g“ d 2 ) 


For small perturbations from the main stream of velocity v 
in the ^-direction, c may be considered equal to the constant 
speed of sound in the undisturbed medium and, in compari- 
son with v, iv = 0j fl*=0, and v z =o. Then 

1 Dp 1 fb 2 <p . o d 2 4> , ,,2 

p Dt~ c 2 V ^ 2 + 


bxbr v bx 1 ) 


With this relation used in equation (A3) the equation for 
the velocity potential may be put in the form given in equa- 
tion (4) of the analysis. 


APPENDIX B 

EVALUATION OF 


s-> 0 


b<f> 


In order to determine the limit of ^ as y-+ 0, it is conven- 


ient to make use of the following substitution: 
2f=(Ti— fi) cos H-fa+fi 


(Bl) 


The expression for 4> (equation ( 10 )) may be written with the 
aid of the following relations (see equation ( 11 )): 


r=p V(r-n) GW) 
= 1 rw ^ g 

P 2 smy 
1 


where 


=0 fo sin 8 

<P=fi P f A(£fi,z+? 0 cos 6) (Ji+f 2 }d8d% (B2) 

Jo Jo 

f fft d/(x-|) , to sin g \ 

c p , Cjg J 
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By the rule for differentiation of a definite integral, 


REFERENCES 


Wy- 


54 

by 


dddi;- f- 


n JT 

, «+JB 


(S3) 


Make use of the following relations: 


bA 


bA_ _ 

by b(z-)-ro cosTj ' 

b^4 bf 0 . 

= cos 8 

bz by 

&(/i+/a) _ cj/i b(t— ti) 

by c)(f— n) dy 


b(g+fo cos 0 ) 

by 


b/ a b(£— t.Q 
"b(f— r 2 ) by 


= ^|(/l~/2)^ sin S 


Then, by integration by parts, the next to the last integral in 
equation (B3) becomes (with cos ddd=dv, ^ ^ (/i+A)— 

+ - fc X? sin " cos »]** 


where the first term vanishes because sin 0=0 at 0=0 and 
0 =tt. Similarly (with sin 8 dd=dv, ^ ^ ^ C/i— /*)=«), 
the last integral in equation (B3) becomes 

«**]>+ 

r«‘ f* ri bA b , , ,, . „ 

»Jo Jo Ue bFbf sm 0 cos 

^ 2 (/ 1 +/ 2 ) cos 2 0j <20<f£ 

where the first term vanishes because/i=/ 2 at 0=0 and 0 =v. 

Then, as y approaches zero from the positive side, there 
results in the limit a contribution only from the first integral 
in equation (B3), 




(B4) 


b£i 


Since ^ changes sign as y changes sign, it follows that 

as y approaches zero from the negative side an equal and 
opposite result is obtained. 
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